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ABSTRACT 

The spectrum of the recently discovered cataclysmic variable star (CV) ASAS-SN 
13cl shows that a secondary star with spectral type K4 (± 2 subclasses) contributes 
roughly half the optical light. The radial velocities of the secondary are modulated on 
an orbital period Porb = 4.86 hr with a velocity semiamplitude K = 246 ± 9 km s“^, 
and the light curve shows ellipsoidal variations and an apparent grazing eclipse. At this 
orbital period, the secondary stars in most CVs are substantially cooler, with spectral 
types near M3. ASAS-SN 13cl therefore joins the small group of CVs with anomalously 
warm secondary stars, which apparently form when the onset of mass transfer occurs 
after the secondary has undergone significant nuclear evolution. 

Subject headings: keywords: stars 


1. Introduction 


Cataclysmic variable stars (CVs) are close binaries consisting of a white dwarf primary that 
accretes matter from a secondary star via Roche lobe overflow. The secondary is more extended 
than the white dwarf, and usually resembles a low-mass main-sequence star. CVs hav e a rich 
pheno menology and have attracted a great deal of observational and theoretical interest. IWarner 
(j 19951 ) gives a dated but useful review of these objects. 


Because the secondary fills its Roche lobe, its mean density is closely con strained by the orbita l 
period, with the secondary’s mean density being larger at shorter periods (jFaulkner et al.lll972l '). 
On the main sequence, the mean density rises toward lower masses, so if the mass-radius relation 
of the secondary is similar to that of main-sequence stars, the orbital period gives a rough proxy 
measurement of the secondary’s mass. Short period CVs therefore tend to have low-mas s secon- 
daries, which tend to have correspondingly cool temperatures and low visible luminosities. iKnigge 


^Based on observations obtained at the MDM Observatory, operated by Dartmouth College, Columbia University, 
Ohio State University, Ohio University, and the University of Michigan. 
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([20061 ) updated the long-known correlation between system’s orbital period and the secondary’s 
spectral type. He found that at a given orbital period, secondaries tended to be near the expected 
main-sequence surface temperature, or somewhat cooler. 


There are exceptions - a few systems are anomalously warm, with K-type secondaries in orbits 
of only a few hours or less. These secondaries appear to have undergone substantial hydrogen 
burning during previous evolution, which increased the mean molecular weight and hence changed 
the mass-luminosit y-temperatur e relat ion. I recently identified an example - CSS 1340 , which 
has Porb = 2.4 hr; iThorstensenI (|2013l ) details this discovery and includes references and further 
background on these objects. 


Here, I report another, less extreme sytem, ASAS-SN 13cl. This previously uncatalogued dwarf 
nova was discovered in outburst by the ASAS-SN survey^on 2013 August 29. Table [U gives further 
details on this object. 


The outburst history of ASAS-SN 13cl is not well-determined. It appears consistently, with 
little obvious variability, on Digitized Sky Survey images served at Space Telescope Science Institute 
1. This ar ea of the sky (at I = 338,6 = —19 deg) is not covered by the Catalina Surveys Data 
Release 2 ( Drake et al. 2009 b 


2. Observations 

All the observations are from MDM Observatory, on Kitt Peak, Arizona. Spectra were obtained 
on 2013 Sept. 11, 13, 14, 15, and 17 UT, and photometric time series were taken on Sept. 16, 17, 
and 19 UT. The system appeared to have returned to quiescence from its recent outburst. 


2.1. Spectroscopy 


The spectra are from the 2.4m Hiltner telescope and modular spectrograph. The instrument 
se tup, observing prot ocols, and reduction procedures were essentially identical to those described 


m 


ThorstensenI ( 20131 ). 


The mean spectrum (Fig. [H top panel) shows the absorption lines of a late-type star together 
with the broad emission lines of hydrogen and Hel commonly seen in dwarf novae at minimum 
light. I estimated the spect ral type and light f r action of the secondary by scaling and subtracting 
spectra of stars classified bv iKeenan fc McNeill (| 19891 ) and looking for the best cancellation of late- 


^The ASAS-SN survey is described at http://www.astronomy.ohio-state.edu/Riassassin/index.shtml 

®The many institutions that made this archive possible are acknowledged at 
http://archive.stsci.edu/dss/acknowledging.html 
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type features. This constrained the spectral type to be K4 ± 2 subclasses, with the secondary 
contributing roughly half the system’s light; the lower trace in the upper panel of Fig. [1] shows the 
result of subtracting a scaled K-star from the observed spectrum. The SDSS magnitud es in Table 
are b roadly consistent with this classification. The i — z color suggests a type near K5 (|Covev et al 
20071 1. with bluer passbands giving earlier types, consistent with an increasing contribution from a 
blue accretion disk toward shorter wavelengths. 

To measure the (absorption) radial velocity of the K star, I cross-correlated the spectra against 
a vel ocity-compensated sum of many late-type spectra, using the rvsao package (|Kurtz &: Mink 
199a). The cross-correlation covered from 5100 to 6500 A, excluding 5800 to 5950 A to avoid Hel 
A5876 and NaD, which can have an interstellar component. Tab le [H gives the velocities. A period 
search of the resulting timeseries with a ‘residual-gram’ method (iThorstensen et al.lll996l ) revealed 
an unambiguous sinusoidal modulation at just under 5 cycle d“^ (Fig. [2]). The Ha velocities, 
measured with a convolution routine, were much noisier than the absorption velocities, but they 
do show significant modulation consistent with the well-defined absorption period. Table [3] gives 
parameters sinusoidal fits to the velocities. The velocity modulation of the absorption lines is 
clearly visible in the lower panel of Fig. [H which shows a phase-averaged greyscale representation 
of the spectra. 


2.2. Time-series Photometry 


The time-series photometry is from an A ndor Ikon DU-937N camera mounted on the MDM 
1.3m McGraw-Hill telescope, as described in iThorstensenI ([2013). Individual exposures were 30 
s, in the V filter. In each frame, the program star, a comparison star, and two check stars were 
measured; the main comparison star was 54 arcsec west and 17 arcsec south of the program star. 
Conditions were generally clear but not photometric, and exposures in which the comparison star’s 
instrumental magnitude more than ~ 0.4 mag fainter than average were excluded from the analysis. 


In the SDSS DRIO (Tabled]), the program object has g = 18.265 and g — r = -1-0.683, while 
the main comparison star has g = 16.194 and g — r = -1-0.545. Using transformations derived by 
Robert Lupton, these correspond to U = 17.866 and 15.875 respectively. Fig. [3| shows the time 
series for the three nights, with the zero-point adjusted to approximate V magnitude, and the lower 
panel of Fig. [2] shows the magnitude plotted against orbital phase. 


The photometry shows variation at twice the orbital period, consistent with ellipsoidal mod¬ 
ulation. Superposed on this is a sharp dip at inferior conjunction of the secondary, which is 
apparently a grazing eclipse. The eclipse is clearly visible on the first and second nights, but only 
just discernible on the third. Table |4| gives the timings of the three eclipses. 
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2.3. Analysis 


Period and Ephemeris. The radial velocities and the grazing eclipses give independent esti¬ 
mates of the period, which agree within their uncertainties. Adopting the weighted mean period 
yields 

BJD eclipse = 2456551.8052(1) 0.20219(8).E, (1) 

where E is an integer. 

Masses. The grazing eclipse requires an orbital inclination i not too far from edge-on, and the 
amplitude of the absorption line velocity variation Aabs should reflect the secondary star’s center- 
of-mass motion with reasonab le hdelity. One is less confident of the dynamical interpretation 


of the Ha velocities (see, e.g., iMarsh et ahl 119871 ): however, the phase at which they cross their 


mean velocity is 0.515 ± 0.025 in the eclipse ephemeris, consistent with the antiphased motion 
expected if they coincide with the white dwarf. If we take the velocities as reliable for dynamics, 
we find q = M 2 /M 1 = K 1 /K 2 = 0.43 ± 0.08. For i = 75 degrees, the mass function then implies 
0.6 Mq < Ml < O. 8 M 0 for the white dwarf, and roughly half that for the secondary. As expecteT 
the secon dary is undermassive for it s spectral type; a recent tabulation by Pecaut and Mamajelo 
(based on IPecant fc Mama)ekll2013l ) estimates 0.7 Mq for typical main-sequence field star at K 6 , 
the late (hence low-mass) end of our spectral type range. 

Distance. The secondary star’s spectral type constrains its surface brightness, and the Roche 
lobe geometry combined with the orbital period lets us estimate the secondary’s radius, with only a 
weak dependence on its mass. The radius and surface brightness yield an absolute magnitude, which 
together with the observed brightness of the secondary (and an estimated interstellar extinction) 
yields a distance. All the quantities entering the calculation are uncertain, and some are correlated 
(e.g., the light fraction attributed to the secondary decreases at cooler spectral types, because the 
metal lines used in estimating the fraction mostly grow stronger toward cooler types). I therefore 
used a Monte Carlo procedure to estimate the distance and its uncertainty, which explicitly included 
the correlation between the secondary’s spectral type and its contribution to the light. For the 
K4 ± 2 secondary, I used magnitudes centered around V = 18.5, and ranging between extremes of 
19.55 (minimum contribution at K 6 ) and 18.10 (maximum contribution at K2). The secondary’s 
mass (which enters weakly) was taken to be uniformly distributed in the range 0.3 ± 0.15. The V- 
ba nd surface br i ghtne ss at each spectral type was computed using a polynomial relation tabulated 
bv iBenermannI (1200611. I ass u med 0.02 < E{B — V) < 0.18, guided by the extinction map of 
Schlegel. Finkbeiner. fc DavisI (j 19981 ). which gives a total extinction of E{B — V) = 0.084 at this 
location. The Monte Carlo procedure yielded a median distance of 1180 pc, with a 68 -percent 
confidence range from 990 to 1400 pc. 

Ellipsoidal Variation. In favorable cases the light variations of the seconda ry star can yield use- _ 

ful constraints on system parameters, especially the inclination. To exploit this. lThorstensen fc Armstrong 


^ available at http://www.pas.rochester.edu/~emamajek/EEM_dwarLUBVIJHK_colors_Teff.txt 
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([2003) developed a code to model the light variation of a Roche-lobe filling star. ASAS-SN 13cl, 
unfortunately, is rather too faint and messy for fruitful light-curve fitting, but nonetheless the 
lower panel of Fig. [2] includes a representative calculatior@ which is useful for comparison. Several 
features stand out: 


1. The sharp dip around phase zero departs signihcantly from the secondary light curve, indi¬ 
cating that it is very likely an eclipse. 

2 . The maxima are unequal, with the second hump (phase 0.75) higher than the hrst. At 
minimum light dwarf novae often show a pre-eclipse b rightening as the ‘hot s pot’ where the 
mass-transfer stream strikes the disk rotates into view (jCoppeians et al.ll2014l present several 
examples of this). While this has the correct sense to explain the inequality, the extra light 
does not persist until eclipse as it does in other systems, so it seems more likely that the 
secondary’s leading and trailing hemispheres differ in mean brightness. 


3. The curve shown has extra light added (assumed constant through the orbit) to approximate 
the contributions from the rest of the system (accretion disk and so on). This tends to flatten 
the light curve and push the system to greater computed distances. The amplitude of the 
computed curve appears to be smaller than that of the data, suggesting that the extra light 
may have been overestimated; however, it is difficult to justify a smaller contribution on the 
basis of the spectral decomposition (Fig. [T]). 


Although the data are noisy and the detailed match imperfect, the modulation’s period and phase 
confirms that it is mostly ellipsoidal. 


Discussion 


The secondary of ASAS-SN 13cl is warmer than expected. At Porb = 4.85 h, the temperature 


and CSS J134052.0-1-151341 ( Thorstensen 


(2nn2a 

), QZ Ser ( 

Thorstensen et al. 

2nn2bb SDSS .1170213.26-^.322954.1 ( 

Littlefair et al. 


200i 




20131 ). but it is signihcant. 


The recent generation of synoptic sky s urveys, most notably the Catalina surveys ([Drake et al 


2009|), ASAS-SN, and the MASTER survey ( Lipunov et al. 201C ) have dramatically accelerated the 




Breedtetal 


2014lb Th e variability criteria favor the 


rate of discovery of new dwarf novae (see 
discovery of large-amplitude outbursts ([Thorstensen fc Skinned l2012l ). In most CVs with periods 


° The curve shown assumes Mi = 0.7 M©, M 2 = 0.35 M 0 , i = 75 deg, Tefi = 4620 K for the M4 secondary 
(following Pecaut &: Mamaiek 2013h . and a contribution from the disk equivalent to an observed V ~ 19.4, which 
includes ~ 0.3 mag of extinction. The model is scaled here to a distance of 1100 pc to match the observational data; 
the distance is nearly identical to that found earlier by an essentially equivalent method. 
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less than a few hours, the secondaries are so faint that the systems can become very faint between 
outbursts, but in systems similar to ASAS-SN 13cl the anomalously bright secondary imposes a 
‘floor’ on the total brightness, which can make the amplitude rather modest. In ASAS-SN 13cl, 
the measured amplitude is less than 3 mag, though it is unlikely to have been caught exactly at 
maximum. At such relatively small amplitudes, discovery is somewhat more difficult. Also, the 
known examples of warm-secondary dwarf novae have relatively long inter-outburst times, or have 
only a single observed outburst. For these reasons, it is likely that a substantial number of these 
warm-secondary systems await discovery. 


I gratefully acknowledge support from NSF grant and AST-1008217. I also thank the MDM 
Observatory staff for their cheerful and excellent support, and the referee for a prompt and useful 
report. 
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Table 1. Properties of ASAS-SN 13cl 


Property 

Value 

Reference^ 

02000 

21*" 38“ 05^046 

SDSS 

<^2000 

+26°38'19”.71 


Outburst date 

2013-08-29.35 

ASAS-SN 

Outburst magnitude 

15.66 


u 

19.097 ±0.028 

SDSS 

g 

18.265 ± 0.008 


r 

17.582 ±0.016 


i 

17.287 ±0.007 


z 

17.055 ±0.014 


V 

17.87 *" 


j 

15.137 ±0.053 

2MASS 

H 

14.814 ± 0.063 


K 

14.701 ± 0.079 



^'References are as follows: SDSS is the Sloan Digital 
Sky Survey Data Rel ease 10; 2MASS is described by 


Skrutskie et alJ (120061 1 . 


^The V magnitude is computed using 

V = g- 0.5784(5 - r) - 0.0038, 


an approximation derived by R. H. Lupton and cited at 

http://www.sdss.org/drT/algorithms/sdssUBVRITransform.html 
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Table 2. Radial Velocities 


Time 

V (absn) 
[km s“^] 

a 

[km s“^] 

V (emn) 
[km s“^] 

a 

[km s“^] 

56546.9338 

-118 

26 

-53 

77 

56546.9408 

-54 

15 

-66 

64 

56548.9009 

-192 

17 

-45 

54 

56548.9093 

-242 

14 

-52 

63 

56549.7790 

-54 

14 

-4 

95 

56549.7875 

22 

11 

-85 

88 

56549.7960 

88 

13 

-123 

no 

56549.8087 

152 

16 

-93 

48 

56549.8172 

206 

15 

-244 

55 

56549.8257 

229 

12 

-185 

52 

56549.8341 

239 

16 

-210 

56 

56549.8426 

226 

11 

-176 

60 

56549.8511 

190 

12 

-240 

59 

56549.8595 

165 

14 

-194 

65 

56549.8680 

152 

10 

-116 

55 

56549.8765 

63 

25 

-132 

52 

56549.8850 

7 

20 

-77 

70 

56549.8935 

-100 

21 

-56 

52 

56549.9020 

-155 

14 

15 

63 

56549.9105 

-144 

16 

-19 

67 

56549.9189 

-220 

14 

33 

92 

56549.9274 

-271 

18 

-65 

97 

56549.9360 

-302 

16 

33 

68 

56550.7524 

-227 

12 

-33 

45 

56550.7609 

-208 

15 

54 

50 

56550.7694 

-176 

13 

-7 

57 

56550.7778 

-152 

17 

-14 

46 

56550.7863 

-102 

13 

-48 

50 

56550.7948 

1 

14 

-139 

61 

56552.8044 

-106 

15 

-4 

57 
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Note. — Times given as barycentric Julian dates of 
mid-integration, minus 2 400 000. The time base is UTC. 
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Table 3. Fits to Radial Velocities 


Description 

To 

(BJD) 

P 

(d) 

K 

(km s“^) 

7 

(km s“^) 

N 

a 

(km s“^) 

Absorption 

56549.7839(11) 

0.20246(20) 

246(9) 

-6(6) 

30 

20 


56549.7841(12) 

[0.20219] 

248(9) 

-7(6) 

30 

22 

Ha emission 

56549.888(5) 

0.2026(9) 

105(17) 

-91(12) 

30 

40 


56549.888(5) 

[0.20219] 

106(17) 

-90(12) 

30 

40 


Note. — Parameters of best-fit sinusoids of the form v{t) = 'y + K sin[27r(t — To)/P]. 
The epoch Tq is expressed as a barycentric Julian day. In the second line for each fit, 
the period is fixed at the weighted mean derived from the radial velocities and eclipse 
timings. 


Table 4. Eclipse Timings 


E 

Barycentric JD 

a 

0 

2456551.8052 

0.0003 

5 

2456552.8164 

0.0007 

15 

2456554.8367 

0.0012 


Note. — E is an integer 
cycle count. Times given are 
barycentric Julian dates of mid¬ 
eclipse. The time base is UTC. 
The last column is the esti¬ 
mated one-sigma uncertainty of 
the timing, in days. 







ASASN-13cl Spectral Decomposition 
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Fig. 1.— Upper panel: Mean spectrum. The individual spectra were shifted to zero velocity 
using the absorption-line ephemeris prior to averaging. The middle trace shows the residual after 
subtracting a scaled spectrum of a K4 star Gliese 570a. The lower trace shows the scaled K4 
spectrum that was subtracted; its zero point is shifted downward by one unit to avoid overlap. 
Lower panel: A greyscale representation the portion of the spectrum from 5200 to 6000 A. Each 
line in the image is derived from a running average of the rectified spectra that lie near the nominal 
orbital phase. The strong absorption line near 5900 A is NaD, and the emission just blueward of 
this is Hel A5876. 
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Fig. 2.— Upper panel: Periodogram of the absorption-line velocities. Middle panel: Absorption 
velocities (blue circles) and Ha (red triangles) emission-line velocities folded on the ephemeris given 
in the text, with the best-fitting sinusoids superposed. Lower panel: V magnitudes, folded on the 
same ephemeris as the middle panel. The solid curve is derived from the ellipsoidal variation model 
described in the text. 
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Fig. 3.— Time series photometry of ASAS-SN 13cl from the McGraw-Hill 1.3m telescope. The 
magnitudes are relative to the main comparison star, with the estimated V magnitude of the 
comparison star (15.875) added to set the zero point, The horizontal bar in the middle panel shows 
the length of the orbital period. 






